Apparatus and method for sampling frequency offset estimation and correction in a wireless communication system

ABSTRACT

A method is described for correcting sampling frequency offset (SFO) of a data packet in a communications system where carrier frequency (f c ) and sampling frequency (f s ) are driven by a common clock source. The method comprises, for each l th  symbol in the data packet, estimating carrier frequency offset (CFO) in a received data packet. From the CFO estimate, an SFO estimate is derived, wherein the SFO is approximately equal to f s  multiplied by said CFO estimate, and divided by f c . An SFO phase correction is generated according to the SFO estimate for each k th  tone in a data stream. The SFO phase correction is then applied to each received data stream.

CROSS-REFERENCE TO RELATED APPLICATION

This application claims the benefit of U.S. Provisional Patent Application 60/691,280, filed Jun. 17, 2005, which is incorporated by reference herein in its entirety.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The invention is generally related to the correction of certain signal impairments in a wireless communication system. More specifically, the invention is directed to the estimation and correction of carrier frequency offset, sampling frequency offset, and phase noise in a wireless communication system.

2. Background

A wireless communication system will typically employ at least one transmitter/receiver pair between which data carrying signals are exchanged. In such a wireless system, the transmitter is typically responsible for encoding the digital information and modulating it onto an analog carrier signal (e.g., a radio signal). Subsequently, the receiver demodulates and decodes the digital information. There are many known techniques for the modulation and demodulation of digital information sent via analog carrier signals.

Regardless of the specific modulation technique used in a given wireless system, the analog carrier signal is subject to certain impairments or degradations during the modulation-transmission-demodulation cycle. Such impairments manifest themselves in a variety of ways including carrier frequency offset (CFO) and sampling frequency offset (SFO). Phase noise is an additional impairment. These impairments, if left uncorrected, can result in corruption of the digital information.

As modulation techniques grow ever more sophisticated, the complexity, speed and sensitivity of wireless communication systems are constantly pushed to the limit. It would be beneficial, then, to continually seek improvements in the estimation and correction of CFO, SFO, phase noise and other such impairments to wireless communication signals.

BRIEF DESCRIPTION OF THE DRAWINGS/FIGURES

The accompanying drawings, which are included to provide a further understanding of the invention, are incorporated in and constitute part of this specification. Together with the description, they serve to provide examples of implementation of the invention. In the drawings:

FIG. 1 illustrates the physical layer of receiver in a wireless communication system.

FIG. 2 more specifically illustrates a portion of the physical layer wherein embodiment of the present invention may operate.

FIG. 3A and 3B illustrate how various signal impairments affect the phase of a signal in a wireless communications system.

FIG. 4A and 4B illustrate a method of estimating and correcting CFO in a received signal.

FIG. 5 illustrates a method of estimating and correcting SFO in a received signal.

The present invention will now be described with reference to the accompanying drawings. In the drawings, like reference numbers indicate identical or functionally similar elements. Additionally, the leftmost digit or digits of a reference number identify the figure in which the reference number first appears.

DETAILED DESCRIPTION OF THE INVENTION

I. Overview of an Exemplary Operating Environment.

While the present invention is described herein with reference to illustrative embodiments for particular applications, it should be understood that the invention is not limited thereto. Those skilled in the art with access to the teachings provided herein will recognize additional modifications, applications, and embodiments that fall within the scope of the claims, and additional fields in which the invention would be of significant utility.

For example, the invention will be described herein in the context of a wireless local area network (WLAN) that uses orthogonal frequency division modulation (OFDM). This exemplary environment was selected because OFDM is a common modulation technique, and is used in the IEEE 802.11a/b/g standards promulgated for wireless local area networks (WLANs). Wireless Ethernet HIPERLAN-2 is also based on OFDM. Additionally, the forthcoming IEEE 802.11n standard for WLAN systems is also expected to employ some form of OFDM. Those skilled in the art will recognize, however, that the methods and systems described herein are equally applicable to other modulation, transmission and demodulation schemes.

An OFDM baseband symbol is the sum of a number of orthogonal sub-carriers, sometimes referred to herein as channels or tones. Each tone is independently modulated by its own data, and thus carries its own data constellation. The modulation schemes may differ, but are typically some type of quadrature amplitude modulation (QAM), complimentary code keying (CCK), or phase shift keying (PSK). Other modulation schemes are within the scope of the invention. A composite baseband signal is typically used to modulate a main radio frequency (RF) carrier having a carrier frequency (f_(c)). OFDM modulation and demodulation are usually implemented with fast Fourier transform (FFT) techniques. When OFDM is used in conjunction with channel coding techniques, it is described as coded OFDM or COFDM.

Each data packet in a WLAN signal using OFDM has a payload portion that includes a number of data symbols. A typical OFDM data symbol uses 20 Mhz of bandwidth and is typically composed of about 52 sub-carriers or tones. There are two types of active tones—data tones and pilot tones. Typically four pilot tones are interspersed among 48 data carrying tones. The pilot tones are typically used as a phase reference. OFDM data packets are typically preceded by a preamble. The preamble can be used for signal detection and synchronization, and typically includes such information as data rate and packet length. Often, the preamble comprises two or more repeated preamble portions. For example, a data packet may be preceded by a short 0.8 μs preamble repeated ten times, followed by a long 3.2 μs preamble repeated two times. The above described data packet and preamble are merely exemplary. Their ultimate form and content are typically determined by whatever standards are implemented in the communications system.

In wide-area wireless broadcasting, one or more receivers can benefit from simultaneously receiving signals from several spatially dispersed transmitters. Such systems may be referred to as a multiple-input multiple-output (MIMO) system, or a multiple-input single-output (MISO) system. MIMO and MISO systems are beneficial because multiple transmitters will destructively interfere with each other only on a limited number of subcarriers or tones, whereas most subcarriers or tones will actually reinforce coverage over a wide area. However, the OFDM signal may nonetheless suffer from time-variations in the channel, which can be manifested as a carrier frequency offset (CFO). Over time, a CFO results in phase rotation of the data symbol. Left uncorrected, the frequency spectrum (and data constellation) will shift, resulting in corrupted data at the receiver output.

Another common impairment is sampling frequency offset (SFO). SFO occurs when a sampling clock is slightly offset (e.g., either slightly faster or slower) from its theoretical or expected speed. For example, a theoretical 80 MHz sampling clock may actually be sampling at 80.08 MHz, or at 79.06 MHz. Such an offset creates problems in data processing. For example, over a thousand samples at 80.08 MHz, the 1000^(th) sample may actually correspond to the 999^(th) piece of data. This is referred to as “sample slip.” However, even over far fewer samples, SFO can effectively cause a small phase roll in the data constellation of a received symbol. Such a phase shift is felt linearly across each received symbol. If left unchecked, such a phase shift can corrupt the received digital information. As explained more fully below, the inventor has determined the relationship between CFO and SFO, thereby permitting effective SFO estimation and correction.

Yet another impairment common to wireless communications systems is phase noise. Phase noise arises as a result of up or down-conversion of the signal by a practical (non-ideal) radio. The effect of phase noise is randomly distributed across a received signal.

As noted above, the inventor has discovered a system and method for improving the estimation and correction of CFO, SFO and phase noise in a wireless communication system. Embodiments of the invention disclosed herein are described in the context of the above described exemplary WLAN environment using an OFDM modulation/demodulation scheme. Estimations and corrections of CFO, SFO, and phase noise are described more fully in the following sections.

II. Carrier Frequency Offset Estimation and Correction

A. Overview

As noted above, a MIMO or MISO wireless communication system employs at least one transmitter/receiver pair between which data carrying signals are exchanged. The transmitter typically encodes the digital information and modulates it onto an analog carrier signal (e.g., a radio signal). In the receiver, the received analog carrier signal, having carrier frequency (f_(c)), is down-converted to baseband, and the digital information carried thereon is demodulated and decoded. The down-conversion, demodulation and decoding takes place in the physical (PHY) layer of a receiver. It is during the demodulation process that any signal impairments must be estimated and corrected to prevent possible corruption of the carried digital information.

CFO estimation and correction can be performed at baseband in the PHY layer of the receiver. The CFO estimation and correction is made using both the preamble and payload portions of the received data packet. For the preamble portion, the correction is iteratively applied to successive preamble portions, and ultimately to the payload portion. For the payload portion, CFO estimates are made and averaged over frequency using multiple tones, and over time using multiple symbols. In MIMO or MISO systems, CFO estimates may also be averaged over space using receive signals from multiple, spatially diverse transmitters. CFO corrections are iteratively applied to the payload as symbols in the data packet proceed through the signal processing pipeline.

CFO estimation and correction is also described in a commonly owned copending application entitled “Apparatus and Method for Carrier Frequency Offset Estimation and Correction in a Wireless Communication System,” attorney docket number 1875.7420001, filed ______, which is hereby incorporated by reference in its entirety.

B. Detailed Description for Carrier Frequency Offset Estimation and Correction.

FIG. 1 is a block diagram of an exemplary PHY layer 100 in a receiver capable of operating in a MIMO or MISO WLAN environment. This environment is exemplary, and was chosen for ease of explanation, and not limitation. PHY layer 100 is coupled to an antenna 105, and receives one or more input analog RF carrier signals 102 a-102 n, upon which digital information has been modulated. PHY layer 100 includes an RF tuner 110, RF filter 115, an analog-to-digital converter 120, an RF demodulator 125, a decoder 130, and a serializer 135. The PHY layer 100 output is coupled to the media access control (MAC) layer. PHY layer 100 also includes a PHY state machine 112, which provides control functions (CONTROL) for certain programmable elements in PHY layer 100. The invention is primarily carried out in a portion of the RF demodulator 125. As noted above, the CFO estimation and correction may be made using both the preamble and payload portions of a data packet. The payload CFO estimate and correction will be described first.

1. Data CFO Estimation

FIG. 2 is a block diagram of a portion of RF demodulator 125. The primary signal path for the received signal begins with a frequency correction block 204, and proceeds through a receive filter and down-sample module 206, a cyclic prefix removal module 208, a time-domain phase locked loop (PLL) correction block 212, a fast Fourier transform (FFT) module 214, a frequency-domain PLL correction block 216, an equalizer 220, and a symbol demapping module 222. For ease of description, the various elements in the primary signal path may be grouped into a first signal processing stage 205, which precedes the FFT module 214, and a second signal processing stage 217, which follows the FFT module 214.

As illustrated, the first signal processing stage 205 precedes the FFT module 214. This stage, the entirety of which is not specifically illustrated, is responsible for down-converting the received data stream to baseband, and preparing the signal for demodulation in the FFT module 214. As shown, the first signal processing stage 205 may include a receiver filter and down-sample module 206, and an optional cyclic prefix removal module 208. As explained more fully below, it is in the first signal processing stage 205 that the preamble portions of the data packet are used to generate a preamble CFO estimate 294. This is accomplished in the preamble estimate module 250. It should be noted that the specific components in the first signal processing stage 205 are exemplary, and are not intended to be limiting beyond what is recited in the appended claims.

The FFT module 214 is primarily responsible for demodulating the received, down-converted data stream. That is, the FFT module 214 expands the received data stream into its sinusoidal components to facilitate further signal analysis, such as symbol demapping and decoding. The time-domain PLL correction block 212 precedes the FFT module 214, while the frequency-domain PLL correction block 216 follows the FFT module 214. As explained more fully below, these PLL correction blocks are typically multipliers that allow frequency or phase corrections to be made to the received datastream as it proceeds through the signal processing portion of PHY layer 100.

The second signal processing stage 217 follows the FFT module 214, and thereafter decodes the impressed digital data for further analysis and distribution in the media access control (MAC) layer of the receiver. As shown, the second signal processing stage 217 may include equalizer 220 and symbol demapping module 222. As with the first signal processing stage 205, the specific components of the second signal processing stage 217 described above are exemplary, and are not intended to be limiting beyond what is recited in the appended claims.

FIG. 2 also illustrates CFO estimator 260. CFO estimator 260 is essentially a phase locked loop (PLL) that receives input from the second processing stages 217, and is coupled in parallel with the FFT module 214. Generally, a PLL is a closed loop frequency control system, that uses feedback to maintain an output signal in a specific phase relationship with a reference signal. The CFO estimator includes a phase detector 262, A pilot/data weighting module 268, a tone average module 263, a symbol-rate loop filter module 264, and a stream average module 266. Each of these modules receives a CONTROL signal from the PHY state machine 112, whereby their operation maybe programmed. As described more fully below, phase detector 262 produces an output signal that is proportional to a phase difference its two input signals. The symbol rate loop filter 264 module may include a voltage controlled oscillator (VCO)/integrator, which is a circuit that produces an ac output signal whose frequency is proportional to the input voltage. The symbol-rate loop filter 264 is programmable circuit used to control the CFO estimator (PLL) dynamics and system performance. One skilled in the art would be familiar with the functions of various PLL circuit elements.

The CFO estimator 260 also includes a number of averaging modules. More specifically, CFO estimator 260 includes a pilot/data weighing module 268, a tone averaging module 263, a and a stream averaging module 266. These programmable averaging modules may be used to further refine the initial CFO estimation generated by phase detector 262.

Generally, CFO estimator 260 accepts three input signals. These inputs, as described more fully below, represent a received signal 292 from second signal processing stage 217, a channel CFO estimate 293 from the channel estimate module 223, and a weighted transmitted signal estimate 296 from multiplexer 226. As explained more fully below, multiplexer 226 assigns an appropriate weight to each transmitted signal estimate (i.e., pilot tone estimate 225 or data tone estimate 240) according to a programmable control signal from the PHY state machine 112. Alternatively, these weights may be assigned by the pilot/data weighing module 268. In an embodiment, pilot tone estimates 225 are weighted more heavily than data tone estimates 240. It should be noted that the particular configuration of the CFO estimator 260, and the manner in which it receives the input signals is merely exemplary. Other configurations could be devised by one of skill in the art with out departing from the spirit and scope of the invention.

To facilitate further description of the operation of CFO estimator 260, the following equations and mathematical conventions are provided. Again, it should be noted that these mathematical conventions are abstract representations, not physical models. As such, they do not represent the actual propagation and impairment of the signal.

The n^(th) transmitter in a MIMO or SIMO system sends a signal impaired by phase noise and CFO as follows: $\begin{matrix} {{x_{n}(t)} = {{\mathbb{e}}^{{j2\pi\Delta}\quad f_{{tx}^{t}}}{\mathbb{e}}^{{j\theta}_{n}{(t)}}{s_{n}(t)}}} & (1) \end{matrix}$ where x_(n)(t) is the received signal, e^(j2πΔf) ^(tx) ^(t)e^(jθ) ^(n) ^((t)) is the phase noise component, and s_(n)(t) is the transmitted signal.

Making the conservative assumption that the phase noise is uncorrelated across a number (N) of transmit antennas, but that the CFO is correlated across the same antennas, the vector representation of the received signal becomes: $\begin{matrix} {\begin{bmatrix} {x_{1}(t)} \\ {x_{2}(t)} \\ \vdots \\ {x_{N}(t)} \end{bmatrix} = {{{\mathbb{e}}^{{j2\pi\Delta}\quad f_{tx}t}\begin{bmatrix} {\mathbb{e}}^{{j\theta}_{1}{(t)}} & 0 & \cdots & 0 \\ 0 & {\mathbb{e}}^{{j\theta}_{2}{(t)}} & \cdots & 0 \\ \vdots & \quad & ⋰ & \vdots \\ 0 & 0 & \cdots & {\mathbb{e}}^{j\quad{\theta_{N}{(t)}}} \end{bmatrix}}\begin{bmatrix} {s_{1}(t)} \\ {s_{2}(t)} \\ \vdots \\ {s_{N}(t)} \end{bmatrix}}} & (2) \end{matrix}$

Imperfect sampling or SFO causes an additional impairment. Assuming that the SFO is correlated across the N transmit antennas, the sampled transmit signal becomes: $\begin{matrix} \begin{matrix} {\begin{bmatrix} {x_{1}\left( {qT}_{s} \right)} \\ {x_{2}\left( {qT}_{s} \right)} \\ \vdots \\ {x_{N}\left( {qT}_{s} \right)} \end{bmatrix} = {{\mathbb{e}}^{{j2\pi\Delta}\quad f_{tx}{qT}_{s}}\begin{bmatrix} {\mathbb{e}}^{{j\theta}_{1}{({qT}_{s})}} & 0 & \cdots & 0 \\ 0 & {\mathbb{e}}^{{j\theta}_{2}{({qT}_{s})}} & \cdots & 0 \\ \vdots & \quad & ⋰ & \vdots \\ 0 & 0 & \cdots & {\mathbb{e}}^{j\quad{\theta_{N}{({qT}_{s})}}} \end{bmatrix}}} \\ {\begin{bmatrix} {s_{1}\left( {q\left( {T_{s} + {\Delta\quad T_{s}}} \right)} \right)} \\ {s_{2}\left( {q\left( {T_{s} + {\Delta\quad T_{s}}} \right)} \right)} \\ \vdots \\ {s_{N}\left( {q\left( {T_{s} + {\Delta\quad T_{s}}} \right)} \right)} \end{bmatrix}} \end{matrix} & (3) \end{matrix}$ where (Ts+ΔTs) is the SFO component.

Over a single symbol, phase noise (e^(jθn(qTs))) can be approximated as constant over one symbols time (e^(jθn,l)). This corresponds to a flat phase shift in the frequency domain. If it is further assumed that the residual CFO for the l_(th) symbol is small, then the CFO can be approximated by a flat phase shift, as illustrated in FIG. 3A. Finally, if it is assumed that the SFO is small, then the SFO can be approximated by a linear phase shift in the frequency domain, as illustrated in FIG. 3B. The frequency-domain representation of the k^(th) tone of the l_(th) transmitted symbol is thus: $\begin{matrix} \begin{matrix} {\begin{bmatrix} X_{1,k}^{l} \\ X_{2,k}^{l} \\ \vdots \\ X_{N,k}^{l} \end{bmatrix} \approx {{\mathbb{e}}^{\frac{{j2\pi\Delta}\quad{f_{tx}{({{2l} - 1})}}T}{2}}{\mathbb{e}}^{\frac{{j2\pi\Delta}\quad{T_{s,{tx}}{({{2l} - 1})}}}{2T_{s}}}}} \\ {\begin{bmatrix} {\mathbb{e}}^{{j\theta}_{1,l}} & 0 & \cdots & 0 \\ 0 & {\mathbb{e}}^{{j\theta}_{2,l}} & \cdots & 0 \\ \vdots & \quad & ⋰ & \vdots \\ 0 & 0 & \cdots & {\mathbb{e}}^{j\quad\theta_{N,l}} \end{bmatrix}\begin{bmatrix} S_{1,k}^{l} \\ S_{2,k}^{l} \\ \vdots \\ S_{N,k}^{l} \end{bmatrix}} \end{matrix} & (4) \end{matrix}$

The receiver phase noise operates on the transmitted signals in a similar way. Thus: $\begin{matrix} \begin{matrix} {\begin{bmatrix} Y_{1,k}^{l} \\ Y_{2,k}^{l} \\ \vdots \\ Y_{N,k}^{l} \end{bmatrix} \approx {{\mathbb{e}}^{\frac{{j2\pi\Delta}\quad{f_{rx}{({{2l} - 1})}}T}{2}}{\mathbb{e}}^{\frac{{j2\pi}\quad k\quad\Delta\quad{T_{s,{rx}}{({{2l} - 1})}}}{2T_{s}}}}} \\ {\left\lbrack \quad\begin{matrix} {\mathbb{e}}^{{j\phi}_{1,l}} & 0 & \cdots & 0 \\ 0 & {\mathbb{e}}^{{j\phi}_{2,l}} & \cdots & 0 \\ \vdots & \quad & ⋰ & \vdots \\ 0 & 0 & \cdots & {\mathbb{e}}^{j\quad\phi_{N,l}} \end{matrix} \right\rbrack} \\ {\left\lbrack \quad\begin{matrix} H_{11,k} & H_{12,k} & \cdots & H_{{1N},k} \\ H_{21,k} & H_{22,k} & \cdots & H_{{2N},k} \\ \vdots & \quad & ⋰ & \vdots \\ H_{{M\quad 1},k} & H_{{M\quad 2},k} & \cdots & H_{{MN},k} \end{matrix}\quad \right\rbrack\begin{bmatrix} X_{1,k}^{l} \\ X_{2,k}^{l} \\ \vdots \\ X_{N,k}^{l} \end{bmatrix}} \end{matrix} & (5) \end{matrix}$

Using the following shorthand notation: $\begin{matrix} {{\mathbb{e}}^{{j\lambda}\quad{k{({{2l} - 1})}}} = {\mathbb{e}}^{\frac{\quad^{{{j2\pi}{({{\Delta\quad T_{s,{tx}}} + {\Delta\quad T_{s,{rx}}}})}}{k{({{2l} - l})}}}}{2T_{s}}}} & (6) \\ {{\mathbb{e}}^{{j\psi}{({{2l} - 1})}} = {\mathbb{e}}^{\frac{\quad^{{{j2\pi}{({{\Delta\quad f_{tx}} + {\Delta\quad f_{rx}}})}}{k{({{2l} - 1})}}T}}{2}}} & (7) \\ {\beta_{n,l} = {\mathbb{e}}^{{j\theta}_{n,l}}} & (8) \\ {\alpha_{m,l} = {\mathbb{e}}^{{j\phi}_{m,l}}} & (9) \\ {\gamma_{{mn},l} = {{\mathbb{e}}^{{j\psi}{({{2l} - 1})}}\alpha_{m,l}\beta_{n,l}}} & (10) \end{matrix}$ the relationship between the impaired received signal vector [Y] and the pure transmitted signal vector [S] (omitting the l subscripts on α, β and γ as well as the k subscript on H_(MN) for readability) is: $\begin{matrix} \begin{matrix} {\begin{bmatrix} Y_{1,k}^{l} \\ Y_{2,k}^{l} \\ \vdots \\ Y_{M,k}^{l} \end{bmatrix} = {{\mathbb{e}}^{{j\psi}{({{2l} - 1})}}{{\mathbb{e}}^{{j\lambda}\quad{k{({{2l} - 1})}}}\begin{bmatrix} \alpha_{1} & 0 & \cdots & 0 \\ 0 & \alpha_{2} & \cdots & 0 \\ \vdots & \quad & ⋰ & \vdots \\ 0 & 0 & \cdots & \alpha_{M} \end{bmatrix}}}} \\ {\begin{bmatrix} H_{11} & H_{12} & \cdots & H_{1N} \\ H_{21} & H_{22} & \cdots & H_{2N} \\ \vdots & \quad & ⋰ & \vdots \\ H_{M\quad 1} & H_{M\quad 2} & \cdots & H_{MN} \end{bmatrix}\begin{bmatrix} \beta_{1} & 0 & \cdots & 0 \\ 0 & \beta_{2} & \cdots & 0 \\ \vdots & \quad & ⋰ & \vdots \\ 0 & 0 & \cdots & \beta_{N} \end{bmatrix}} \\ {\begin{bmatrix} S_{1,k}^{l} \\ S_{2,k}^{l} \\ \vdots \\ S_{N,k}^{l} \end{bmatrix}} \end{matrix} & (11) \\ \begin{matrix} {= {{\mathbb{e}}^{{j\psi}{({{2l} - 1})}}{{\mathbb{e}}^{{j\lambda}\quad{k{({{2l} - 1})}}}\begin{bmatrix} {\alpha_{1}\beta_{1}H_{11}} & {\alpha_{1}\beta_{2}H_{12}} & \cdots & {\alpha_{1}\beta_{N}H_{1N}} \\ {\alpha_{1}\beta_{2}H_{21}} & {\alpha_{1}\beta_{2}H_{22}} & \cdots & {\alpha_{2}\beta_{N}H_{1N}} \\ \vdots & \quad & ⋰ & \vdots \\ {\alpha_{M}\beta_{1}H_{M\quad 1}} & {\alpha_{M}\beta_{2}H_{M\quad 2}} & \cdots & {\alpha_{M}\beta_{N}H_{MN}} \end{bmatrix}}}} \\ {\begin{bmatrix} S_{1,k}^{l} \\ S_{2,k}^{l} \\ \vdots \\ S_{N,k}^{l} \end{bmatrix}} \end{matrix} & (12) \end{matrix}$

In terms of phase offset (γ_(mn)), the following relationship emerges for a MIMO system: $\begin{matrix} {\begin{bmatrix} Y_{1,k}^{l} \\ Y_{2,k}^{l} \\ \vdots \\ Y_{M,k}^{l} \end{bmatrix} = {{{\mathbb{e}}^{{j\lambda}\quad{k{({{2l} - 1})}}}\begin{bmatrix} {\gamma_{11}H_{11}} & {\gamma_{12}H_{12}} & \cdots & {\gamma_{1N}H_{1N}} \\ {\gamma_{21}H_{21}} & {\gamma_{22}H_{22}} & \cdots & {\gamma_{2N}H_{2N}} \\ \vdots & \quad & ⋰ & \vdots \\ {\gamma_{M\quad 1}H_{M\quad 1}} & {\gamma_{M\quad 2}H_{M\quad 2}} & \cdots & {\gamma_{MN}H_{MN}} \end{bmatrix}}\begin{bmatrix} S_{1,k}^{l} \\ S_{2,k}^{l} \\ \vdots \\ S_{N,k}^{l} \end{bmatrix}}} & (13) \end{matrix}$ where the [Y] vector comprises received signals 292 at the M^(th) receiver, for the k^(th) tone, and for the l_(th) symbol. The [H] vector represents an estimated signal impairment (e.g., CFO or SFO) between a receiver/transmitter pair (M,N), with γ representing the phase offset. The [S] vector comprises the transmitted signal estimate 296 from the Nth transmitter, for the l_(th) symbol and, for the k^(th) tone.

If a vector [Z] is established as being the product of [H] and [S] from equation 13, and we solve equation for phase offset γ, the phase offset γ may be represented by the inverse [Z] matrix multiplied by vector [Y]. γ≈1/[Z]·[Y]  (14)

It should be noted that the above mathematical description provides a convenient model of CFO, SFO and phase noise. However, they are not physical models, and therefore do not represent the actual propagation and impairment of the signal. Nonetheless, the above described mathematical model is useful for illustrative purposes and will be referred to at points in the following description.

With equation 14 in mind, a general method of estimating and correcting CFO in a receiver operating in the above described exemplary WLAN environment is now described in FIGS. 4A and 4B. As noted above, the received data packet includes multiple data streams. Each data packet has a preamble having a number of preamble portions. The preamble is followed by payload that includes a number of data symbols L. Each l_(th) data symbol, in turn, has a number of tones K. Each k^(th) tone carries its own unique data. In a MIMO or SIMO system, there are a number M of received data streams in each data packet.

Generally, the methods 400A and 400B involve, for each data symbol, making both a preamble CFO estimate and a per-tone CFO estimate. An initial CFO correction is made based on the preamble CFO estimate. Then, for an l^(th) symbol in a data packet, the per-tone CFO estimate is averaged across frequency, time and space to generate a data CFO estimate, P_(data). CFO in the receiver is then updated using the P_(data) estimate, and/or the preamble CFO estimate, P_(preamble). The steps are repeated for each l^(th) data symbol in the data packet.

More specifically, method 400A begins with step 402, where a preamble CFO estimate, P_(preamble) _(—) _(stream(m)), is made by correlating repeated preamble portions that precede each data packet. The derivation of a preamble estimate is described more fully below in section 2. In step 403 a, a weight is assigned to each P_(preamble) _(—) _(stream(m)). Thereafter, a single preamble CFO estimate, P_(preamble), is generated by averaging the M weighted per-stream P_(preamble) _(—) _(stream(m)) CFO estimates, according to step 403 b.

An initial CFO correction is then made based on the single preamble CFO estimate P_(preamble), according to step 404. Next, according to step 406, a per-tone CFO estimate is successively averaged over frequency, time and space to derive a data CFO estimate P_(data). Step 406, as explained more fully below, is repeated for each l^(th) data symbol. Finally, according to step 408, any CFO in each received data stream is then corrected according to the data CFO estimate, P_(data), and/or the preamble CFO estimate P_(preamble).

FIG. 4B describes in more detail how the CFO estimate is successively averaged over frequency, time and space to generate a final CFO estimate that is applied to the receiver. For ease of description, method 400B is described with reference to the mathematical conventions described above, in view of FIG. 2. Nonetheless, as noted above, these mathematical conventions and reference to FIG. 2 should not be construed as limiting beyond what is recited in the appended claims.

The method 400B begins with step 410 where, for an l^(th) symbol in a data packet, a per-tone CFO estimate P_(tone)(k) is generated for each of the K data tones in the data symbol. This per-tone CFO estimate is based on (i) a received signal [Y]; (ii) a channel estimate [H]; and (iii) a data estimate or pilot information [S]. As illustrated in FIG. 2, a control signal (CONTROL) to an exemplary multiplexer 226 provides one method of selecting whether the data 240 or pilot 225 information is to be combined with the channel estimate [H] 291 and the received signal [Y] 292. Use of the multiplexer 226, however, is merely exemplary, and one of skill in the art could envision other methods for selecting between the data and pilot information, as necessary.

A weight is selectively assigned to each of the K per-tone CFO data estimates in accordance with step 415. The tone weights are based on channel conditions. As noted above, individual tones may be data tones, or pilot tones. Data tones carry the digital data constellation that has been modulated onto them, while the pilot tones carry predetermined, known data constellations. Thus, the estimates for the transmitted tones may be made from either data tones or pilot tones.

If the estimate is made from the data tone, then it is taken from the output of symbol demapping module 222, which uses well known decision directed techniques to estimate the contents of the data tone. Such known techniques include making a “nearest neighbor” estimate on individual data points in a data constellation, or by providing “decision boundaries” for each point in the data constellation. If the tone estimate [S] is made from the pilot tones, then the estimate is taken from a pilot generator 224, which simply provides the known data for the pilot tone.

In an embodiment, the pilot tone estimates are weighted more heavily than the data tone estimates. In an exemplary embodiment, multiplexer 226 assigns an appropriate weight to each transmitted signal estimate according to a programmable control signal CONT from the PHY state machine 112. However, one of skill in the art could envision other methods for assigning an appropriate weight to the K per-tone CFO estimates. For example, a pilot/data weighting module 268 could be provided after phase detector 262 to assign appropriate weights to the per-tone CFO estimates.

In an embodiment, the weighted transmitted tone estimate 296 [S] is received into phase detector 262, along with a channel CFO estimate 293 [H]. These signals may be combined in phase detector 262—for instance, with a multiplier—and may be represented mathematically by vector [Z] in equation 14. As such, this vector would represent a refined, per-tone, transmitted signal estimate. Also received at phase detector 262 is received data signal 292. Received data signal 292 may be represented by vector [Y] in equation 14, and can be taken from the output of frequency domain PLL correction block 216, prior to the estimates made by in symbol demapping module 222.

As illustrated in FIG. 3A, the phase offset due to CFO is essentially constant across individual data symbols. Thus, according to step 420, a per-symbol CFO estimate, P_(sym)(l) is next generated for the l^(th) symbol by averaging the K weighted per-tone CFO estimates. This is accomplished by tone average module 263. The per-symbol CFO estimate represents an average CFO over the frequency bandwidth occupied by the K tones in the data symbol.

In steps 425 and 430, a second weighted average is calculated. According to step 425, a weight is selectively assigned to each of the l per-symbol CFO estimates generated in step 420. Next, according to step 430, the weighted l per-symbol CFO estimates are averaged to generate a data stream CFO estimate, p_(stream)(m). The data stream CFO estimate, p_(stream)(m), thus represents an average CFO over time as data packets proceed through the signal processing pipeline. The averaging is accomplished by symbol-rate loop filter module 264. In an embodiment, the most recent per-symbol CFO estimates are successively given the highest weight, while the weight given to older per-symbol estimates are successively diminished, each time the data stream CFO estimate, p_(stream)(m), is calculated.

If there are multiple received data streams M, as in a MIMO or SIMO communications system, then steps 410-430 are repeated for each M data stream to generate M data stream CFO estimates, p_(stream)(m), according to steps 435 and 440. In step 435, a weight is assigned to each m per-stream CFO estimate. The relative weights are assigned based on channel conditions using, for example, the signal quality in the channel. A higher weight would be assigned to channels with a higher signal quality. According to step 440, the M data stream CFO estimates are then averaged to generate a single data CFO estimate, P_(data). The data CFO estimate, P_(data), essentially amounts to averaging the CFO over space—i.e., over two or more spatially diverse transmitters. The averaging is accomplished by stream average module 266.

Recall from above that the data packet has both a data portion and a preamble portion. The data CFO estimate, P_(data), generated above in steps 410 to 440 may be used with the preamble CFO estimate to generate or update a packet CFO estimate, P_(packet). Depending on channel conditions and system behavior, the packet CFO estimate, P_(packet), may comprise either the preamble CFO estimate, the data CFO estimate, or a combination of both. Ideally, the most reliable CFO estimate—typically the preamble CFO estimate—would be used. However, a noisy environment or a bad channel could dictate using the data CFO estimate, or a combination of both. Thus, according to step 445, a packet CFO estimate, P_(packet), is generated using the preamble CFO estimate, and/or the data CFO estimate, P_(data).

Finally, according to step 450, the CFO in each data stream is corrected according to the packet CFO estimate P_(packet). The correction is applied by the CFO correction generator 280, which receives the output of CFO PLL 260 and applies it to the time domain PLL correction block 212. Correction block 212 can be implemented as a multiplier that, in an embodiment, precedes FFT module 214. As can be seen, although the CFO estimate is accomplished with demodulated data (i.e., post FFT data), the frequency correction is applied prior the FFT operation. Steps 410 through 450 are repeated for each data symbol in the data packet, according to step 455.

In sum, for an l^(th) data symbol in a data packet, a per-tone CFO estimate is made for each of the K tones in the data symbol, and successively averaged over frequency, time and space to generate a data CFO estimate. The data CFO estimate may be combined with a preamble CFO estimate (see below) to update a packet CFO estimate. The packet CFO estimate is then applied to the data stream passing through the signal processing pipeline. In this fashion, after only a few data symbols have passed through the signal processing pipeline, any carrier frequency should be substantially eliminated, thereby permitting accurate demapping and decoding of the digital data impressed thereon.

2. Preamble CFO Estimation and Correction

As part of the overall system CFO estimation and correction process, a preamble CFO estimate may be made using the data packet preamble. Each data packet is preceded by a preamble. The preamble has repeated preamble portions. For example, a data packet may be preceded by a short 0.8 μs preamble repeated ten times, followed by a long 3.2 μs preamble repeated two times. The ultimate form and content of the preamble are typically determined by whatever standards are implemented in the communications system. Regardless of the particular form of the preamble, the contents of the repeated preambles are substantially identical. This allows repeated preambles to be compared to each other over time and, if a phase roll is detected, a rough estimate of CFO may be obtained.

To make a preamble CFO estimate, repeated preamble portions are correlated. For ease of explanation, the correlation may be described mathematically as follows. Assume that a first received data signal X₁ may be represented by X₁=S₁e^(j2π(fc)(Ts))+V₁. Where S₁ is the transmitted sample, f_(c) is the combined transmitted and received center frequency of the analog carrier signal, T_(s) is the sampling rate, and V₁ is a noise term. A second received sample of a substantially identical portion of a successive repeated preamble portion may be represented by X_(N+1)=S₁e^(j2π(fc)(N+1)(Ts))+V₂. A correlation of the samples yields X₁*X_(N+1)=|S₁|²e^(j2π(fc)(Ts)N)+V₁·S₁e^(j2π(fc)(N+l)(Ts). As can be seen, the phase argument in each of the correlated terms is 2π(f) _(c))(Ts)N. Thus, if the correlation of successive samples is summed over repeated portions of the preamble, that summation will be essentially equal to the phase argument 2π(f_(c))(Ts)N. Because N and Ts are known quantities, one can solve for actual f_(c) and compare it to the theoretical f_(c) to make a preamble estimate of CFO.

In a MIMO or SIMO system, there are multiple, substantially identical received data streams from spatially diverse transmitters. This allows the preambles in the M data streams to be correlated and summed in the same fashion as repeated preambles in a single data stream. Furthermore, the preamble based estimates for each of the M data streams may be weighted according to received signal quality estimates before combining into a single CFO estimate. The preamble CFO estimates from received signals having a better signal quality would be weighted higher than preamble CFO estimates from received signals having a lower signal quality. This allows further refinement of the preamble CFO estimate in a MIMO or SIMO system.

The preamble CFO estimate can be applied to successive portions of the data packet. The successive portion of the data packet may be the payload portion. Alternatively, the successive portion of the data packet may be a successive preamble portion. Either way, such a correction permits progressively narrowing CFO in the network as the data packets proceed through the first signal processing stage 205.

III. Sampling Frequency Offset Estimation and Correction.

As described above, another common impairment is sampling frequency offset (SFO). SFO occurs when a sampling clock is slightly offset (e.g., either slightly fast or slow) from the desired or theoretical sampling rate. For example, a theoretical 80 MHz sampling clock may actually be sampling at 80.08 MHz, or at 79.06 MHz. Such an offset creates problems in data processing. For example, over a thousand samples at 80.08 MHz, the 1000^(th) sample may actually correspond to the 999^(th) piece of data. This is referred to as “sample slip.” Sample slip may occur in either a positive or negative direction, depending on the cumulative direction of phase roll. However, even over far fewer samples, SFO can effectively cause a small phase roll in the data constellation of a received symbol. As illustrated in FIG. 3B, such a phase shift is felt linearly across each received symbol. If left unchecked, such a phase shift can corrupt the received digital information, especially at the outer ends of the symbol bandwidth.

The inventors recognized a relationship between CFO and SFO that permits effective SFO estimation and correction. In a communications system such as the exemplary OFDM WLAN described above, there is typically a single crystal oscillator clock from which other system clocks are derived.

Therefore, in such a system, the same crystal oscillator drives both the sampling frequency (f_(s)) and the carrier frequency (f_(c)). In such systems, the SFO can be expected to be of a substantially identical proportion of the theoretical f_(s), as CFO is of the theoretical f_(c). Put differently, if CFO and SFO were to be expressed in parts-per-million, the their value would be substantially identical. This relationship permits effective estimation and correction of SFO, where CFO is known.

Returning to FIG. 2, a system for estimating and correcting SFO is illustrated. As fully described above, FIG. 2 shows a preamble estimate module 250 that uses the preamble portion of the data packet to make a preamble CFO estimate. FIG. 2 also shows the CFO estimator 260, which further refines the CFO estimates by, for each symbol, first making a per-tone CFO estimate, and successively averaging it over time, frequency and space (if necessary) to further refine and produce a data CFO estimate. The output of the CFO estimator 260 thus represents a data CFO estimate, P_(data). The above described method of estimating CFO may be implemented in the present invention. However, the present invention is not so limited, and other methods of estimating CFO known in the art may be used as well.

FIG. 2 further illustrates an SFO correction module 270 that is coupled to the output CFO estimator 260. SFO correction module 270 includes an SFO calculator 272 that uses the assumed relationship between CFO and SFO to generate an SFO estimate. As described above, SFO is approximately equal to the theoretical f_(s) multiplied by the CFO estimate 261 (from CFO estimator 260) and divided by the theoretical f_(c). An SFO phase correction generator 274 is coupled to SFO calculator 272. SFO correction generator 274 receives the SFO estimate, generates an SFO phase correction, and applies it to the received data stream. In an embodiment, the SFO phase correction is applied to the frequency domain PLL correction block 216, which follows FFT module 214, to achieve a phase roll in the opposite direction as the estimated SFO. In an alternate embodiment (not illustrated), the SFO phase correction may be applied to the data stream as a time correction (ΔT_(s)) by using a linear interpolator. Such a correction would typically be applied prior to FFT module 214. A time shift by ΔT_(s), ensures that the data is sampled at ideal points on the data constellation.

SFO phase correction generator 274 also includes an SFO rotor 276. SFO rotor 276 is used to track overall phase roll. As noted above, if overall phase roll due to SFO reaches a certain threshold value, sometimes referred to as a “guard interval,” sample slip can occur. In other words, as the SFO causes the linear phase roll across symbols, the data or pilot tones at the edges of the symbol bandwidth become offset to a degree that is beyond the decision boundary of symbol demapping module 222. At a certain threshold value (typically approximately one half the sample spacing (T_(s)/2)) the SFO phase roll tracker (PRT) 276 sends a signal 277 to cyclic prefix removal module 208, instructing it to either advance or retard a data sample. This may be accomplished by rotating the data signal by 360 degrees in the desired direction. For example, the sample may be advanced by rotating the phase of the received data stream by negative 360 degrees when the phase roll is greater, in a negative direction, than the predetermined guard interval. Alternatively, the sample may be retarded by rotating the phase of the received data stream by positive 360 degrees when the phase roll is greater, in a positive direction, than the predetermined guard interval. As illustrated in FIG. 2, cyclic prefix removal module 208 follows the receive filter and downsample module 206.

FIG. 5 illustrates a method 500 for estimating and correcting SFO in a received data stream. As with the CFO correction method, method 500 is accomplished for each l^(th) symbol in a data packet. As noted, the described method presumes that f_(c) and f_(s) are driven by a common clock source, such as a crystal oscillator. According to step 510, an estimate is made of the CFO in the received data packet. Two methods of estimating CFO are fully described above. However, other methods of estimating CFO are known in the art, and may be used as well.

In step 515, an SFO estimate is derived from the CFO estimate. As described above, where f_(c) and f_(s) are driven by a common clock source, the SFO can be expected to be of a substantially identical proportion of the theoretical f_(s), as CFO is of the theoretical f_(c). More specifically, SFO is approximately equal to the theoretical f_(s) multiplied by said CFO estimate, and divided by the theoretical f_(c). The SFO estimate may be derived from either the preamble CFO estimate 294, the data CFO estimate 261, or another known method for estimating CFO.

An SFO phase correction is generated in step 520 for each k^(th) tone, according to the derived SFO estimate. Given the above described relationship between CFO and SFO, a determination is made on how far, and in which direction, a corrective phase roll is to be applied to the data streams passing through the primary digital signal processing pipeline.

In an embodiment, and according to step 525, SFO is corrected in each received data stream using the frequency domain PLL correction block 216 that follows FFT module 214. Application of the SFO phase correction causes a phase roll in each received data stream that corrects small SFOs. Alternatively, an interpolation filter (not shown) could be employed to shift the sample in time according to the derived SFO estimate. This would have the effect of ensuring that the sample points match up to ideal points in the data constellations embodied in the plurality of data tones. Such an interpolation filter would typically be employed prior to FFT module 214. In any event, the SFO phase correction generator 274 ensures that the SFO correction is in an appropriate form, depending on the type of correction is applied.

According to step 530, phase roll due to SFO is tracked. As described above, phase roll due to SFO grows over time. If phase roll due to SFO becomes too large, sample slip may occur. Therefore, if phase role due to SFO passes a certain predetermined threshold value, then the phase of the received data stream will be rotated 360 degrees in the appropriate direction to prevent sample slip. As noted above, the threshold value is approximately equal to one half the sample interval T_(s).

Thus, according to step 535, a determination is made as to whether the phase roll due to SFO is greater than the predetermined threshold. If the phase roll is not greater than the threshold, then sample slip is not imminent, and SFO PRT 276 continues to monitor phase roll, according to step 530. However, if the phase roll does exceed the predetermined threshold, then, according to step 540, the data stream is advanced or retarded by one sample, as appropriate. Simultaneously, the SFO phase correction is rotated by 360 degrees in the appropriate direction. If the sample is advanced, then the SFO phase correction is rotated negative 360 degrees. If the sample is retarded. Then the SFO phase correction is rotated positive 360 degrees. After the sample has been either advanced or retarded, the system returns to tracking phase roll due to SFO according to step 530.

In sum, SFO can be estimated from CFO where both f_(c) and f_(s) are driven by a common clock source. An SFO correction is derived from the SFO estimate, and applied to the data stream. Phase rotation due to SFO is also tracked. If the phase rotation due to SFO surpasses a threshold, then the sample is either retarded or advanced to prevent SFO induced sample slip. The SFO phase correction is simultaneously rotated by 360 degrees, as appropriate, and the system continues to track phase roll due to SFO.

Various embodiments of the present invention have been described above, which are capable of being implemented on an interactive graphics machine. It should be understood that these embodiments have been presented by way of example only, and not limitation. It will be understood by those skilled in the relevant art that various changes in form and details of the embodiments described above may be made without departing from the spirit and scope of the present invention as defined in the claims. Thus, the breadth and scope of the present invention should not be limited by any of the above-described exemplary embodiments, but should be defined only in accordance with the following claims and their equivalents. 

1. A method for correcting sampling frequency offset (SFO) of a data packet with one or more data streams (M) in a communications system where carrier frequency (f_(c)) and sampling frequency (f_(s)) are driven by a common clock source, the method comprising: for each l^(th) symbol in the data packet: (a) estimating carrier frequency offset (CFO) in a received data packet; (b) deriving a SFO estimate from said CFO estimate, wherein the SFO is approximately equal to f_(s) multiplied by said CFO estimate, and divided by f_(c); (c) generating a SFO phase correction according to said SFO estimate for each k^(th) tone in the data packet; and (d) applying said SFO phase correction to said one or more data streams.
 2. The method of claim 1, wherein said received data packet has a preamble portion, and wherein said CFO estimate is derived from said preamble portion.
 3. The method of claim 1, wherein said received data packet has a payload portion, and wherein said CFO estimate is derived from said payload portion.
 4. The method of claim 1, further comprising (e) tracking phase roll due to SFO; (f) when said phase roll is greater, in a negative direction, than a predetermined threshold, advancing a start position of the sample processing by one sample thereby preventing sample slip, and simultaneously rotating the SFO phase correction negative 360 degrees; and (g) returning to step (e).
 5. The method of claim 1, further comprising (e) tracking phase roll due to SFO; and (f) when said phase roll is greater, in a positive direction, than a predetermined threshold, retarding a start position of the sample processing by one sample thereby preventing sample slip, and simultaneously rotating the SFO phase correction positive 360 degrees; and (g) returning to step (e).
 6. The method of claim 1, wherein said received data packet is demodulated using a fast Fourier transform (FFT).
 7. The method of claim 6, wherein said SFO correction is a phase correction that is applied after said FFT.
 8. The method of claim 6, wherein said SFO correction is a time correction that is applied before said FFT.
 9. The method of claim 1, wherein said CFO estimate is generated, for an l^(th) data symbol in the data packet, by deriving a per-tone CFO estimate and successively averaging said per-tone CFO estimate over frequency, time and space to derive the CFO estimate. 